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Abstract 

Building on the technological success of the IceCube neutrino telescope, we outline a prospective low-energy extension 

that utilizes the clear ice of the South Pole. Aiming at a 10 Mton effective volume and a 10 MeV threshold, the detector 

would provide sufficient sensitivity to detect neutrino bursts from core-collapse supernovae (SNe) in nearby galaxies. 

The detector geometry and required density of instrumentation are discussed along with the requirements to control 

the various sources of background. We find that the resulting detector will be able to detect SNe from beyond 10 Mpc, 

fSj delivering between 11 and 46 regular core-collapse SN detections per decade. It would further allow to study more 

^ speculative phenomena, such as optically dark (failed) SNe, where the collapse proceeds directly to a black hole, at a 

Qh detection rate similar to the regular SNe. We find that the biggest technological challenge lies in the required large 

^^ number of large area photo-sensors, with simultaneous strict limits on the allowed noise rates. If both can be realized, 

^— ,. the detector concept we present will reach the required sensitivity in a cost effective manner and hence offers a route to 

future routine observations of SNe with neutrinos. 

^H Keywords: supernovae, core-collapse, neutrino physics, neutrino detectors 
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In 1987, Supernova 1987A (SN 1987A) exploded in the 
Large Magellanic Cloud at a distance of only 50kpc, lead- 
ing to the first detection of neutrinos from outside our solar 
system. Despite the fact that only ~20 supernova neutri- 
nos were detected in total pL^, a wealth of papers has been 
published in its wake (see e.g. [2] for a summary), reflecting 
the numerous and fundamental roles that neutrinos play in 
astrophysics as well as in particle physics (e.g. see [3l, l4l). 

Given today's detectors, a supernova in our Galaxy 
would result in ~ 10^ neutrino events detected individu- 
ally in Super-Kamiokande [5 and other large low-energy 
neutrino detectors, as well as up to millions of neutrinos 
detected through an increase in noise rate in IceCube [6 . 
However, with an expected rate of only 1-3 Galactic SNe 
per century [71 [8] , the chance for a detection during the 
lifetime of the experiments is not overwhelmingly large. 
In the fortunate case of a SN detection, the uniqueness of 
the progenitor system will make it difficult to disentangle 
the astrophysical diversity from the effects due to parti- 
cle physics (e.g. neutrino oscillations) that will impact the 
light curve and energy spectra. 

As pointed out in [7 and [8 , the situation will change 
drastically once neutrino detectors reach the sensitivity 
threshold to detect "mini-bursts" of neutrinos from su- 
pernovae in neighboring galaxies. Not bound to our own 
Galaxy, the rate of SN observations will depend only on 
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the size of the detectors. As we will show in section |5| an 
effective volume of ~ 10 Mtons is sufficient to detect SNe 
at a rate of ~ 1—4 per year — albeit most of them with 
less than ten individual neutrino events. Despite the low 
number of detected neutrinos, these routine observations 
would provide a wealth of information and allow entirely 
new studies [8]. What follows is a brief and incomplete 
summary of the scientific benefits of a large supernova neu- 
trino detector. 

The total SN rate in our local universe would be de- 
termined in a novel and less biased way. Given an appar- 
ent mismatch in rates — only about half the expected rate 
of SNe is actually detected by optical surveys [9 — a di- 
rect measurement would hence allow to solve the riddle of 
missing supernovae. 

Furthermore, a significantly sensitive detector allows 
to test models predicting additional neutrino bursts, such 
as failed or dark (i.e. optically unobservable) SNe [10], 
merger of binary neutron stars [11] or the formation of 
quark stars [12 . Dark SNe are core-collapse objects that 
directly form a black hole (BH) . Electromagnetic radiation 
is strongly suppressed, since the photons don't have time 
to escape and are swallowed by the forming BH. Neutri- 
nos, on the other hand, can escape, and the expected burst 
from such an event is both more luminous and hotter [10^ . 
Average neutrino energies can be roughly twice as large as 
in the case of the collapse to a neutron star (NS) and hence 
open the opportunity to identify the collapse to a BH. An- 
other opportunity is the identification of questionable SN 
candidates, e.g. luminous blue variables [13 , as "super- 
nova impostors" by the non-detection of neutrinos [8 . 
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In addition, detections of neutrino bursts can be used 
to trigger early optical or X-ray observations. Having a 
precise timing for the moment of explosion and observ- 
ing the shock-breakout in electromagnetic radiation, will 
allow to infer a wealth of information about the progeni- 
tor system. In the absence of a direction for the SN, the 
follow-up could focus on observation of nearest galaxies, 
since these are the only ones for which a neutrino detec- 
tion is expected. 

With a larger number of supernova neutrinos, even a 
broader physics program can be accessed. Due to the com- 
plexity of the involved processes, modeling of supernova 
explosions is still a challenge today and has significant 
variance in e.g. the predicted mean neutrino energy [M]. 
Determination of the neutrino luminosity and energy spec- 
trum will provide valuable input to these models. 

From the observation of the arrival times of the neu- 
trinos from SN 1987A, a limit of tti^ < 5.8 eV has been 
set on the mass of the electron neutrino [T5] . By observ- 
ing neutrinos from supernovae with higher statistics, more 
stringent limits could be set. In addition, by its impact on 
the predicted flux, the mass hierarchy can be addressed 
as well, provided a sufficient number of neutrinos is ob- 
served [6j. 

Motivated by this scientific potential, several megaton 
scale neutrino detectors are currently planned (e.g. Deep- 
TITAND [8 , Hyper-Kamiokande [16, IX and UNO [18,). 
Those are either water Cherenkov detectors located in 
mines or marine detectors, similar to ANTARES. In this 
paper, we explore the potential of a ~ 10 Mton detector 
in the Antarctic ice shield. In the existing IceCube detec- 
tor that has been optimized for TeV-PeV energies [19 , SN 
neutrino bursts are typically searched for by looking for a 
collective enhancement of photomultiplier noise rates [6] . 
Due to the large sensor spacing and consequently high en- 
ergy threshold, attempting to detect individual neutrino 
events in this configuration significantly reduces the effec- 
tive mass and hence distance at which SNe can be de- 
tected [20 . A dedicated effort is now under way to reduce 
the energy threshold to a few GeV in the PINGU (Phased 
IceCube Next-Generation Upgrade) project [21 , with one 
of its goals being the determination of the neutrino mass 
hierarchy through the MSW effect [22]. Building on this 
effort and the expertise accumulated with IceCube and 
AMANDA, we explore the capabilities of a dedicated low- 
energy extension to study individual supernova neutrino 
events. The challenge is to reduce the energy threshold of 
the experiment by three orders of magnitude while control- 
ling the background at a level required for the detection 
of supernovae. 

While we will focus on this aspect in this paper, it 
should be noted that such a detector will not be limited 
to the detection of neutrinos at a few MeV, but will also 
provide unprecedented sensitivity at the GeV-scale. Not 
only the neutrino oscillation sector will be accessible with 
increased precision — it will also provide sensitivity to 
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Figure 1: Time evolution of v^^ flux for a SN at 1 Mpc 
in the Lawrence-Livermore model [26 (full line represents 
the differential flux and dashed line the integral). 



other astrophysical phenomena such as collisional heating 
in gamma-ray bursts |23l |24] . 

The paper is organized as follows: The SN models that 
we use to benchmark a future detector are discussed in sec- 
tion |2J followed by a description of our detector simulation 
and the optimization of the detector configuration in sec- 
tion [3J Section [4] is a discussion of the various background 
sources that can be anticipated. Finally, we give an esti- 
mate on the expected SN rate in section |5] and conclude 
with a discussion of the results in section |6l 



2. Neutrinos from core-collapse supernovae 

In core-collapse supernovae, neutrinos are produced in 
course of the formation of a neutron star which follows the 
gravitational collapse of a massive progenitor star. Elec- 
trons and protons rapidly fuse to neutrons, emitting elec- 
tron neutrinos, which results in a sharp (~ 10 ms) delep- 
tonization peak in the neutrino flux which marks the birth 
of the neutron star. However, most supernova neutrinos 
(~ 90%) originate from thermal pair creation [25]. These 
neutrinos, released over a characteristic diffusion time of 
1-10 s (cf. Fig.jl]), are of all flavors, with flavor ratios, light 
curve and energy spectrum still being under debate [M]. 

An early simulation of neutrino production in SN ex- 
plosions has been provided in the so-called Lawrence Liv- 
ermore (LL) model [26]. In the LL model, the neutrino 
spectrum is parametrized by 
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where L is the luminosity and {E) the average energy while 
P determines the width of the spectrum. The parameters 
for the LL simulation as well as for a competing model by 
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Table 1: Parameters of the Lawrence Livermore (LL) 
model [2 , the Thompson, Burrows, Pinto (TBP) 
model ^27j and the dark supernova model [10 for the neu- 
trino spectrum of a core-collapse supernova. The table 
presents the stellar mass of the progenitor Mq , the aver- 
age neutrino energy {E^J, the pinch parameter ^^^ (see 
Eqn.fTl) and the time-integrated luminosity in anti-electron 
neutrinos Lj^^. The TBP model does not lead to an ex- 
plosion, i.e. no luminosity emerges from the simulation. 
Instead, the luminosity from the LL model is assumed. 
Please note that for dark SNe Eqn. [l]is not valid, instead 
the positron spectrum given in ^ was used. 



Thompson, Burrows and Pinto (TBP) [27] are reproduced 
in Table [H While the latter model does not result in an 
explosion of the supernova, [14 suggests that the resulting 
luminosity should be of the same order of magnitude. A 
large variety of other neutrino emission models exists [28^ , 
but we restrict ourselves to these two commonly used ones 
for ease of comparison. Alternative models predicting low- 
energy neutrinos include dark supernova models in which 
stars heavier than ~ 25 M© form black holes. It is assumed 
that if the black hole rotates fast enough, it explodes as 
a hypernova, otherwise the event should be very faint or 
even dark in optical emission, while even more luminous 
in neutrino flux than ordinary supernovae [29l |30]. The 
resulting neutrinos have average energies of (Ej^) ~ 20 — 
24 MeV [10 . 

The neutrinos from the different collapse scenarios can 
best be detected via inverse beta decay 
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that requires a threshold energy of Eiy > 1.806 MeV. We 
use the approximation for the energy dependence of the 
cross-section called "Naive +", which is presented in [3T1. 
Since the cross-section rises with energy, the resulting positron 
spectrum is harder than the initial neutrino spectrum. 
Figure [2] shows the resulting positron spectra for the two 
different supernova models as well as for a model of dark 
supernovae. Neutrino oscillations, which harden the spec- 
trum, are taken into account only for the dark model, for 
which we directly take the positron spectrum from [10]. 

For a typical positron energy of 20 MeV the corre- 
sponding track length is ~ 10 cm, resulting in ~ 3600 
Cherenkov photons [6 . Since the light yield scales lin- 
early with the positron track length, and hence with the 
positron energy, the average amount of light produced per 
neutrino is model dependent. 



Figure 2: Positron spectrum for the Lawrence-Livermore 
(LL) [26 , Thompson, Burrows, Pinto (TBP) [27 and a 
dark supernova model [10 for a SN at 1 Mpc and an effec- 
tive detector volume of 1 Mton. 



3. Detector simulation and optimization 

In this section we describe the simulation and opti- 
mization of a possible detector in the Antarctic ice capa- 
ble of detecting MeV neutrinos with high statistics. Be- 
ing embedded in natural ice, special attention has to be 
given to the optical properties of the detection medium 
and hence the location of the detector. Two parameters 
affect the propagation of the Cherenkov light - the absorp- 
tion length Xa and the effective scattering length Ag. Both 
parameters have been subject to dedicated measurement 
campaigns during the installation of the AMANDA and 
IceCube detectors [32H34] . From these measurements, two 
natural choices for a location of the detector arise: 

1. 750 - 1050 m where the AMANDA- A [35^ detector is 
located. We will call this diffusive ice. 

2. 2150 — 2450 m where the DeepCore array is located. 
We will call this clear ice. 

In the diffusive ice in a depth range around 900 m, 
the absorption length of up to A^ ~ 350 m is exception- 
ally large. However, the presence of air bubbles results 
in very strong scattering with a scattering length of only 
Ag ~ 0.3 m [36^. This results in an effective propagation 

length of Xp = y^^^ ~ 6m after which the photon flux 

has dropped by a factor e~^ [32 . While the photons can- 



not travel large distances, they are confined within a small 
volume for a rather long time of A^/c ~ 1 /is before finally 
being absorbed. This leads to a large detection proba- 
bility, as long as the light is emitted in the vicinity of a 
photo-sensor. In the diffusive ice the light propagation is 
well described analytically as a random walk with absorp- 
tion [32 , allowing for convenient simulation of the photon 
tracking. 



Much larger scattering lengths of Ag ~ 20 — 50 m are 
found in the clear ice, at depth ~ 2300 m, where air bub- 
bles have fully degenerated due to the high ambient pres- 
sure and only a small dust concentration was measured. 
The absorption length in this clear ice is A^ ^ 20 — 90 m, 
depending on wavelength [33", 34 , providing a good envi- 
ronment for a neutrino detector. In particular, the large 
scattering length allows for precise reconstruction of the 
event position and direction. However, compared to diffu- 
sive ice, the photons are more likely to travel through the 
ice without hitting a photo-sensor. For the clear ice, the 
Photonics code [ 37H39] — a Monte Carlo simulation devel- 
oped for IceCube — is employed for accurate simulation of 
the photon propagation. It includes the full depth depen- 
dence of the scattering and absorption properties of the 
ice [33] . 

In what follows, detectors in both depths will be dis- 
cussed. We will show that in a pure photon counting argu- 
ment the diffusive ice is superior to the clear ice. However, 
the possibility to reconstruct individual events in the clear 
ice will allow to address a broader spectrum of physics 
questions. Also in terms of background, the diffusive ice 
is disfavored. In the clear ice, the atmospheric muon flux 
is greatly suppressed due to the large overburden and our 
studies show that reconstruction of neutrino direction - 
impossible in the diffusive ice - will be important for dis- 
crimination of solar neutrinos. We will also show that re- 
construction of the neutrino events is required to overcome 
the large background from self- noise in the photo- sensors. 
Hence, the clear ice is the preferred location. 

In these first simulations, each of the detectors is as- 
sumed to consist of 61 strings of 300 m length, with one 
photo-sensor per meter. In order to get a fairly homoge- 
nous distribution of sensors in the ice, the strings are ar- 
ranged in a filled hexagon, similar to the layout of the Ice- 
Cube detector. For the diffusive ice, we model the photo- 
sensors as vertically aligned tubes of 100% efficiency with 
/ = 1 m length and tt ■ d = 1 cm circumference. Averaging 
over all directions and considering project ional effects, the 
cylindrical sensor has an effective area of (Aeff) = j-l-d = 
25 cm^. For the simulation of the clear ice, we assume 
IceCube-like photo-sensors, i.e. the so-called digital opti- 
cal module (DOM) [40 . Each DOM includes a 10" Hama- 
matsu PMT, which is integrated into a pressure resistant 
glass sphere that also includes the electronics for HV gen- 
eration and in-ice digitalization of the PMT signal. Ice- 
Cube employs both regular efficiency and high-quantum 
efficiency (HQE) PMTs. The dark noise rate of individ- 
ual DOMs averages around 500 Hz [6j. IceCube DOMs 
have a non-trivial directional sensitivity [41] which is in- 
corporated in the Photonics simulation package [39] and 
hence included in our simulation. From laboratory mea- 
surements [42 , we obtain a Cherenkov- weighted effective 
area of 19.4 (26.3) cm^ for an IceCube DOM equipped with 
a regular (HQE) PMT. In order to achieve a comparable 
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Figure 3: The effective detector mass at trigger level for a 
61-string geometry (averaged over the LL model spectrum) 
as function of string spacing. In clear ice (dashed) the 
effective mass is several times lower than in diffusive ice 
(solid) due to ice properties. Only if the sensor effective 
area is increased by a factor of four (dotted) , similar results 
can be obtained. 



detector sensitivity in the clear and diffusive ice, we find 
that we need to boost the photo-sensitive area for the de- 
tectors in clear ice. We therefore scale the sensor response 
for the DOMs we simulate in the clear ice to obtain ef- 
fective areas of 35.6 cm^ and 4 x 35.6 cm^ per meter of 
instrumented string. 

In order to optimize the detector geometry, the effec- 
tive neutrino detection volume is calculated as function of 
the lateral spacing between strings. We employ a Monte 
Carlo method where neutrino vertices are generated ho- 
mogeneously within and beyond the detector volume. The 
positron energy is sampled according to the respective en- 
ergy spectrum shown in Fig. [2J The average number of de- 
tected photons is calculated for each sensor with respect to 
the neutrino vertex. For the diffusive ice, this is done via 
analytic functions, for the clear ice, the Photonics code 
is employed. The actual number of detected photons is 
drawn from a Poisson distribution and the arrival times 
are sampled for each of these photons from the respective 
arrival time distribution. 

We then calculate the effective mass^ i.e. the mass 
of the geometrical simulation volume multiplied with the 
fraction of neutrinos that are detected above threshold. 
This threshold is set at a minimum of five photo-sensors 
being hit by photons. The threshold of five sensors is 
chosen to allow reconstruction of the vertex position and 
positron direction which add five degrees of freedom. Fig- 
ure [3] shows the resulting effective mass as a function of the 
string spacing for the clear and diffusive ice. Generally one 
finds that, as the string distance increases, the geometri- 



cal detector volume increases as well, but the fraction of 
detected neutrino events with at least five hit sensors de- 
creases, yielding a maximum effective volume for a string 
spacing around 30 m for a Lawrence-Livermore neutrino 
spectrum. The maximum effective mass at trigger level in 
the diffusive ice is ^ 9.5Mton. In the clear ice, a much 
smaller effective mass of ~ l.SMton is achieved using the 
same effective photo-sensitive sensor area as in the diffu- 
sive ice. In order to achieve the goal of a close to 10 Mton 
detector, we increase the sensor area by a factor four to 
arrive at a maximized effective mass of ~ 8 Mton in the 
clear ice. 

We have also considered an alternative cylindrical de- 
tector layout similar to the Super-Kamiokande [5] or Hyper- 
Kamiokande [161 E] layouts by simulating a detector that 
consists of only one ring of strings. The strings form a 
cylindrical wall that encompasses an un-instrumented vol- 
ume. However, after optimizing the ring radius for the 
same number of sensors and strings, we generally find an 
effective volume that is reduced by 45% compared to the 
filled hexagonal geometry. This is mainly because the 
cylinder is open at the top and bottom so that many 
events are lost. Deploying a horizontal layer of sensor in 
the antarctic ice is not possible in a cost-efficient manner. 
Hence, we do not explore the cylindrical geometry any fur- 
ther. 

4. Background studies 

Up to this point, we have solely implied the trigger 
requirement — that the number of sensors hit by photons 
be at least five for each neutrino event. In this section, we 
will discuss the dominant sources of background and how 
they can be controlled by implying additional constraints 
on the distribution of photon hits. Contributions from 
random noise, atmospheric neutrinos, atmospheric muons 
and solar neutrinos are considered. 

Neutrinos from SNe come in bursts. To be distin- 
guished from uncorrelated background or noise triggers, 
a SN will need to produce a certain multiplicity of neu- 
trino triggers within a given time window. To estimate 
the allowed background trigger rate, we assume that a su- 
pernova detection can be claimed if three or more events 
occur within a time window of AtsN ~ 1 — 10 s. The rate 
of background triggers that we can allow in order to expect 
not more than one fake supernova claim per year is 
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3.99 mHz for A^^at = 1 s 
0.86 mHz for A^^at = 10 s 



(3) 



Hence, a background trigger rate of at most a few mHz 
is allowed if one wants to effectively detect SNe via the 
outburst of neutrinos. 



4.1. Sensor noise 

Since the Antarctic ice shield is a very low-radioactivity 
environment, the main sources of random noise are in- 
troduced by the detector itself: radioactive isotopes and 
thermal noise in the photo-sensors and electronics. For 
the IceCube modules, this results in a dark count rate of 
~ 500 Hz [6 . While it is known that some fraction of 
the noise is not purely random, but correlated in time, we 
will neglect this for simplicity of the discussion presented 
here. As was shown above, the effective photosensitive 
area will have to be increased to about 0.014 m^ per in- 
strumented meter to achieve the effective target mass of 
10 Mtons. While conventional PMTs of this size are avail- 
able, the required thickness of a glass pressure vessel of this 
size and the cost to drill holes of the corresponding diame- 
ters make this option very unattractive. New photo-sensor 
technologies (e.g. based on wavelength shifters as light col- 
lectors) are currently discussed for deployment in future 
extensions to IceCube, that offer to increase the effective 
photo-sensitive area while simultaneously significantly re- 
ducing the noise rate. These technologies are still in the 
design phase, so we use dark noise values of 500 Hz, 50 Hz 
and 10 Hz as templates in absence of a measurement. 

As shown above (Eqn. |3|, the rate /noise of fake neu- 
trino triggers caused by random noise hits has to be below 
^ 1 mHz if only one fake supernova event is tolerated per 
year. We will now calculate this noise trigger rate depend- 
ing on the total number of modules in the detector TVtot, 
the random noise rate per module /^n, the number of hit 
modules ntrig that is required for a neutrino event to trig- 
ger and the trigger time window ttrig. Assuming that one 
module has registered a random noise hit and opened the 
trigger window, the probability Pm for any module to also 
see at least one noise hit during this time window ttrig is 
complementary to the probability to register no hit, which 
can be calculated according to a Poisson distribution 



Prr. = 1 



(4) 



The probability for a noise event Pnoise is the probability 
that at least ntrig — 1 more modules also encounter a noise 
hit in the time window. It can be calculated with the bino- 
mial distribution, again via the complementary probability 
of registering less than ntrig — 2 noise hits: 
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where Bcum{m \ n,p) = YT=o il) P^ (f -p)"""^ is the cu- 
mulative binomial probability for up to m successes out of 
n tries with probability p. Apart from boundary effects, 
the rate of noise triggers in the detector is then 



/n 



■ fm ■ Nu 



(6) 



Using generic values for the dark count fm = 500 Hz, 
A/tot = 18300 for a detector with 61 strings and 300 mod- 
ules per string and requiring ntrig = 5 hits in ttrig = 
1000 ns the rate of fake SN events is /noise = 9 MHz, i.e. 
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Figure 4: Dimensions and effective sensor density of the 
final detector layout in the clear ice with 127 strings 
and 25 m string spacing as compared to the IceCube and 
PINGU/DeepCore. 



is 9 orders of magnitude above the allowed value from 
Eqn. [3| Even assuming a module dark noise rate as low 
as fm = 10 Hz, still /noise = 7.4 Hz, well in excess of what 
can be tolerated. This clearly shows that it is necessary 
to apply intelligent trigger algorithms that take advantage 
of the non-uniform distribution of photons from neutrino 
interaction and thus limit the number of modules consid- 
ered by the trigger. 

For the diffusive ice, this will be basically impossible as 
the scattering will smear out the photon distribution while 
at the same time the long photon survival time requires 
trigger time windows of ttrig > 2 /is to catch a major frac- 
tion of signal hits. Hence, low fake SN event rates in the 
diffusive ice can only be obtained at the cost of extremely 
low trigger efficiencies. For the clear ice in contrast, the 
event geometry (i.e. the distribution of light throughout 
the detector) is well preserved. In the following, we present 
two strategies that reject fake events induced by sensor 
noise. As these lead to reduction of the detection effi- 
ciency, we have extended the detector to a configuration 
with 127 strings in the same hexagonal geometry to retain 
the lOMton target set initially. This final geometry and 
the corresponding effective sensor density (in units of Ice- 
Cube regular efficiency modules) are illustrated in Fig. [4] 
and will be used throughout the following. 
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Figure 5: Distribution of signal hits in a simulated de- 
tector: spatial distance Ar to the vertex position vs. time 
difference At between the hit and the neutrino interaction. 



4.1.1. Phase space cut 

In case of a neutrino event, a positron from inverse 
beta decay produces Cherenkov light along its track, which 
is short compared to the scattering length in ice. These 
events are thus characterized by a localized emission of 
light. In contrast, events generated by random noise have 
no intrinsic topology and have hits distributed uniformly 
throughout the detector. The rate of noise events can thus 
be significantly reduced by only considering the part of the 
detector where a signal event would create most of the hits. 

For this we define a fiducial volume within the detec- 
tor that depends on the time and position of the neutrino 
vertex. As that vertex is not known a priori, we use a 
X^-minimization of the time residual of the photon prop- 
agation time minus the expected photon arrival time for 
straight propagation from the vertex to reconstruct the 
vertex position. Using this simple method, a positional 
resolution of '^ 15 m can be achieved. To incorporate this 
limited knowledge of the vertex position, a random Gaus- 
sian smearing with 10 m standard deviation in each coor- 
dinate is applied to the vertex position of our simulated 
events. Fig. [5] shows the distribution of hits in the two- 
dimensional phase space given by the time between the 
hit and the neutrino interaction and the distance between 
the hit module and the reconstructed vertex. The 10% 
and 90% quantiles of the distance are calculated for this 
phase space (green lines in Fig. pj red line is the median) , 
giving the inner and outer radius of a spherical shell ex- 
panding with time and containing ^ 80% of the signal hits. 

We now reassess the noise trigger rate by limiting the 
hits that contribute to an event to those within this phase 
space region (cf. Fig. [6| for all times up to tmax- As the 
phase space volume (i.e. the number of sensors that can 
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Figure 6: Illustration of the phase space cut and the lo- 
cal coincidence cleaning on a noise event. Each colored 
sphere corresponds to a hit optical module. The blue re- 
gion indicates the expanding spherical shell used in the 
phase space cut, dashed circles indicate the radius of the 
local coincidence cleaning. 
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contribute to the trigger) quickly increases with tmax, 
fake SN event rate /noise is a strong function of tmax as 
well. We choose the value of tmax so that it limits /noise ^ 
1 mHz for a given sensor noise rate /^ , number of hits 
required to form a trigger ntrig, and the string spacing 
c^str of the detector and calculate the effective mass for SN 
signal neutrinos for each configuration. 

We simulate the detector in the clear ice as described 
above with 127 strings of 300 m length (c.f. Fig.[4|. Fig.JT] 
shows the resulting effective mass after this phase space 
cut is applied and Tab. [2] lists the optimal parameters as 
function of module noise /^. While this method reduces 
the trigger efficiency for neutrino interaction by up to a 
factor of 10, at the same time the fake trigger rate is re- 
duced by many orders of magnitude. For the module dark 
noise rate of 500 Hz provided by IceCube DOMs, achiev- 
ing a fake trigger rate of 1 mHz is only possible at the cost 
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Table 2: Optimal parameters for string spacing ((igtr), re- 
quired hit multiplicity (ntrig) and maximum hit time win- 
dow (tmax) together with resulting effective mass Meff as 
function of module noise rate jm- The parameters are 
chosen so that the phase space cut limits the rate of noise 
triggers to /^Jg^ = 1 mHz. The noise rates for the same 
configuration without the phase space cut are shown for 
comparison. 



of simultaneously increasing the number of hits required 
for a trigger, decreasing the string spacing and using a 
short trigger time window. Due to these extreme choices, 
the trigger efficiency for neutrino events is low and the 
effective mass of this detector configuration falls short of 
the initial target of 10 Mton. In order to retain 10 Mton, 
reducing the sensor self-noise rate to ~ 10 Hz per meter 
instrumented string is essential. 

^.1.2. Local coincidence cleaning 

While taking advantage of the localization of the event, 
the above method does not account for the fact that light 
from the positron is emitted in a preferred direction on the 
Cherenkov cone. As demonstrated for IceCube, requiring 
a local coincidence between photon hits is a very efficient 
way to reduce the effect of random noise [43^. A hit is 
required to be accompanied by at least another hit within a 
certain radius Ar and time window At in order to fulfill the 
local coincidence criterion. Fig. |6] shows an illustration of 
this process. We have applied such a radius-time cleaning 
to the simulated events and find it to provide additional 
background rejection power on top of the phase space cut. 
In order to optimize the parameters At and Ar, we apply 
the hit cleaning method to signal and noise events and 
then re-optimize the parameters of the phase space cut for 
the cleaned data set. The reduction of noise hits allows to 
increase tmax for the same /noise of 1 mHz. Fig.[8]shows the 
effective mass obtainable with this additional cleaning for 
10 Hz module noise to be increased by 60% from 7.5 Mton 
to 12 Mton. 

J^.l.S. Conclusion 

Operation of a very dense detector at the lowest pos- 
sible trigger threshold introduces a strong dependence of 
the trigger rate on the module noise. Algorithms exploit- 
ing the spatial and temporal distribution of the light (that 
can not be used in the diffusive ice) have to be employed 
to obtain a sufficiently low rate of fake events. While even 
more advanced algorithms considering the full event topol- 
ogy would still allow a moderate increase of the effective 
mass, this study demonstrates the importance of low-noise 
sensors. Correlated noise that was not treated here will be 
more difficult to reject and provides additional motivation 
to seek low dark noise rates in future photo-sensor R&D. 

4.2. Muon background 

Muons crossing the detector or passing nearby are eas- 
ily separable from the SN neutrino signal via the huge 
amount of Cherenkov light produced by the extended track. 
An additional outer layer of the detector, naturally pro- 
vided by IceCube in the deep clear ice, would ensure that 
even muons passing by at large distances or stopping just 
above the detector can be identified as such. 

For a conservative first estimate of the effect of atmo- 
spheric muons, we assume all muons reaching the top of 
the detector being energetic enough to cross it. From [44] 
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Figure 7: Effective mass of the simulated 127-string detec- 
tor at trigger level (solid) and after a phase-space cut is 
applied to limit the fake trigger rate to 1 mHz for module 
noise rates of 500 Hz (dashed), 50 Hz (dashed, circles) and 
10 Hz (solid, circles). 



Figure 8: Effective mass of the simulated 127 strings de- 
tector with 10 Hz noise rate per module at trigger level 
(solid), with phase space cut applied (dashed, circles) and 
with phase space cut re-optimized after the radius-time 
cut is applied (solid, circles). 



we obtain a muon flux of $^ = 8 • 10~^ cm~^s~^sr~^ for 
the top of the clear detector at 2150 m depth, giving a 
muon passing rate of about 



Ru = ^u • 300 m X 300 m • TT ?^ 230 Hz 



(7) 



using a detector cross-sectional area of (300 m)^ and an 
effective solid angle of n that accounts for the lower flux 
from angles closer to horizon where the ice shield is thicker. 
A muon traveling through the entire detector has a track 
length of ~ 300 m and emits about A^o ^ 10^ Cherenkov 
photons on its path {^ 360 photons per cm). We assume 
as a worst case that all photons are trapped within the 
detector volume by scattering to compute how long the 
photons will remain detectable within the detector before 
it is absorbed. The number of photons after a path x = 
ct/n is given by: 



N^ = No- e 



(8) 



with tabs = "^^a/c = 1.5 • 10~^ s, for an absorption length 
of \a ^ 100 m (as found in the clear ice) and a refractive 
index of n = 1.31. After a time 



t = iabs m — « 7 /Lis 



(9) 



the number of unabsorbed photons in the detector due to 
a crossing muon has fallen below A^^ = 1 which is well 
below our trigger threshold. We therefore conservatively 
assume each passing muon to illuminate the detector for a 
time interval of 7 /is, which has to be accounted for as dead 
time for supernova neutrino detection, corresponding to a 
fraction of R^ • 7/is = 0.16% of detector operation time. A 
detector in the diffusive ice will not only experience a much 



higher muon passing rate of over 5000 Hz as it is located 
much closer to the surface. Because of the long absorption 
length of ~ 350 m, photons will remain in the detector for 
~ 25 /is, causing an effective dead time of ~ 14% during 
which the detector is saturated with photons. 

These passing muons are easily identified because of the 
dense instrumentation and can be rejected applying a veto 
on down- going particle tracks. However, muons also pro- 
duce spallation products by fragmenting ^^O nuclei [45] . 
The decay of these products mimics low-energy neutrino 
events and is a serious background. Super-Kamiokande 
has demonstrated [45l [46] that a cut on a likelihood func- 
tion including spatial and temporal distance from the pass- 
ing muon as well as energy loss of the muon can be used 
to efficiently remove spallation events. However this re- 
sults in an additional 20% effective dead time for Super- 
Kamiokande which is located at a depth of 2700 meters 
water equivalent [45] , comparable to the deep location con- 
sidered here. For the shallow diffusive ice around 900 me- 
ters depth in contrast, spallation backgrounds may render 
neutrino detection entirely impossible. 

4.3. Solar neutrino background 

Put aside the cosmic neutrino background — which is 
too low in energy to be detectable — the by far dominant 
flux of neutrinos at the earth stems from the sun, where 
neutrinos are abundantly produced in several different fu- 
sion cycles [2 . As only u^ and no Uq are generated in the 
sun, these solar neutrinos cannot undergo inverse beta- 
decay. The dominant interaction for solar neutrinos is elas- 
tic scattering on electrons. Interactions on oxygen nuclei 
are about two orders of magnitude less important [4T and 
are ignored here. Furthermore, only ^B and hep neutrinos 



need to be considered, since all other solar neutrino fluxes 
are below the detector's energy threshold at ~ 1.7 Me V. 
Since the flux of hep neutrinos is three orders of magni- 
tude lower than the ^B flux they are neglected here as 
well [13 [49]. 

To calculate the interaction rate, we use an analyti- 
cal expression for neutrino-electron elastic scattering from 
Eqn. (5.25) in [2] 
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with Tg being the kinetic energy of the recoil electron in 
the laboratory frame and 



ctq 
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gi = - -\- sin Ow — 0.73 92 = sin Ow — 0.23 

and fold it with the energy-dependent effective mass of 
our detector given in Fig. |9jb). Taking the shape of the 
^B neutrino spectrum from [50 and normalizing it to a 



total flux of ^{^B) 



2.44 • lO^cm-^s 



[47], 



we arrive 



at an approximate solar neutrino event rate of ~ 75 mHz 
for a detector in the clear ice at trigger level. Additional 



application of the phase space cut (cf. Sec 4.1.1) and lo- 
cal coincidence cleaning (cf. Sec 4.1.2) reduces this rate 
to ~ 28 mHz, still significantly larger than the maximum 
allowed rate of random background events /^^ = 1 mHz. 
However, a number of methods can be applied to fur- 
ther discriminate the rate of solar neutrino events. The 
bulk of solar neutrinos is less energetic than the bulk of 
SN neutrinos (cf. Fig. [91(a)). Changing the trigger require- 
ment from 5 hits to 7 hits increases the energy threshold 
and thus reduces the expected event rate of solar neutrinos 
by a factor of nearly four, while reducing the signal effi- 
ciency for SN neutrinos from the LL model by only 37%. 
Furthermore, the electron emerging from the elastic scat- 
ter roughly keeps the direction of the incident neutrino, 
while the inverse beta decay effectively randomizes the 
positron direction. For a sufficiently dense instrumented 
array such as Super-Kamiokande, an angular resolution 
of - 30° at Ee = lOMeV is feasible [461 [47]. Assuming 
a one-sigma (68% of the events) angular window of this 
size to reject neutrinos from the direction of the sun, we 
estimate to be able to further discriminate the solar neu- 
trino rate to /© = (1 - 0.68)^ • 28 mHz ^ 1 mHz, while 
retaining (1 - (30°)V^s/c2y)^ > 93% of the SN events. The 
most powerful rejection will be achieved using a likelihood 
method that incorporates both the direction and energy 
for each event within the time window. In absence of a 
full reconstruction we refrain from a more detailed discus- 
sion, but note that the large rate of solar neutrinos can be 
significantly reduced. It is also worth mentioning that a 
fundamental prerequisite is the ability to reconstruct the 
direction of the emerging lepton, which can only be done in 
the clear ice where most photons are only slightly deflected 



5 hits 


0.753 


0.0038 


3.760 


5 hits, PS+RT 


0.282 


0.0032 


2.513 


6 hits, PS+RT 


0.158 


0.0030 


2.032 


7 hits, PS+RT 


0.079 


0.0027 


1.586 



Table 3: Average number of events per 10 s from solar 
and atmospheric neutrino background (up to 100 MeV) 
as well as the SN neutrino signal for a LL supernova in 
10 Mpc distance for different numbers of hit modules and 
with noise cleaning (PS+RT) applied. 



before absorption. In the diffusive ice, all directional in- 
formation is lost during the random walk of the photons. 

4.4- Atmospheric neutrino background 

Cosmic rays colliding with the Earth's atmosphere pro- 
duce Ue and Ue in similar abundance. The dominant com- 
ponent are the electron anti-neutrinos, interacting via in- 
verse beta-decay with a cross-section two orders of mag- 
nitude higher than the neutrino-electron elastic scattering 
of electron neutrinos [46 . Taking the atmospheric neu- 
trino flux calculations from 51 , the neutrino-electron 



elastic scattering cross-section from Eqn. 10 and the cross- 
section for inverse beta-decay given in the phenomenolog- 
ical parametrization in [31J ("Naive +" model) and in- 
tegrating the event rate from 3 — 100 MeV, we arrive at 
an expected trigger level event rate of 0.004 mHz for u^ 
and 0.4 mHz for u^ triggering on 5 hit sensors in the clear 
ice. The resulting spectrum, also shown in Figure [91(d), 
peaks well above the peak of SN interactions, allowing to 
further discriminate these events. While somewhat more 
abundant, elastic scattering of i^^ is only possible via neu- 
tral current interactions, and thus has a factor six smaller 
cross-section. 

Another component of the background are invisible 
muons that are produced by low-energy atmospheric muon 
neutrinos. These are muons below the threshold for Cherenkov 
light emission, that decay to an electron that has - due to 
its lower mass - a velocity above the Cherenkov thresh- 
old. These Michel electrons have been measured by Super- 
Kamiokande [46j _52^ and amount to ~ 90 events per year 
in their effective volume of 22.5ktons. Scaling this up to 
^ 4000 events per Mton per year, we arrive at a corre- 
sponding rate of ~ 1.5 mHz for a 12 Mton detector. While 
already small compared to the background of solar neutri- 
nos, further reduction can be achieved by using the sur- 
rounding IceCube detector to veto accompanying atmo- 
spheric muons in the deep ice. 

4.5. Summary 

The large sensor multiplicity requires intelligent trig- 
ger and selection algorithms to cope with the backgrounds 
arising from sensor self-noise, atmospheric muons, solar as 
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Figure 9: (a) The flux for solar [4T, "50] and atmospheric [5T] neutrinos as function of energy. The Supernova neutrino 
flux on Earth according to the LL model [26 , normalized to an average power of 4.9 • 10^^ ergs~^ in v^^ for a supernova 
distance of 10 Mpc is shown assuming a burst duration of 10 s. (b) Effective mass as function of lepton energy for the 
127-string detector in clear ice, triggering on 5 hit modules, after applying noise cleaning. Above 100 MeV, a linear 
extrapolation (dashed) is used, (c) The effective area for a 127-string detector for elastic scattering and inverse beta- 
decay, (d) Event rate per MeV as function of neutrino energy for solar, atmospheric and supernova neutrino fluxes from 
(a) with effective area from (c). 
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Figure 10: SN detection probability for the simulated 127 
strings detector after application of noise cleaning as de- 



scribed in section 4.1 Shown are results for the LL model, 



the TBP model, and a dark supernova model (compare Ta- 
ble fl]), for the detection of at least either 3 or 10 neutrino 
events from the SN. 



well as atmospheric neutrinos. In this challenging back- 
ground situation, a detector in the diffusive ice close to 
the surface is strongly disfavored due to a) the inability to 
reject self-noise triggers by their hit distribution pattern, 
b) the much higher muon background leading to significant 
downtime, c) the inability to use the surrounding IceCube 
and DeepCore detectors as muon veto and d) the solar 
neutrino backgrounds which cannot be vetoed effectively 
due to the lack of directional reconstruction. Even in clear 
ice, the suppression of noise to a sufficiently low level will 
be a major challenge and contingent on future improve- 
ments in sensor development. Due to the larger scattering 
length compared to the diffusive ice, the hit pattern is 
less localized, but contains directional information, mak- 
ing the application of event reconstruction both necessary 
and possible. 

5. Expected supernova detection rate 

Knowing the detector effective area as a function of 
neutrino energy (Figure ^(c)), we can proceed to calculate 
the sensitivity to a supernova neutrino burst with a given 
spectrum. The probability to detect a supernova is calcu- 
lated from Poisson statistics. We consider a supernova as 
detected if at least three neutrino events trigger the de- 
tector within 10 s, which is the threshold we can expect if 
all other backgrounds can be controlled to within 1 mHz. 



Figure [To] shows the SN detection probability as function 
of distance to the SN, using the three considered models, 
for the 127-string detector in the clear ice described in sec- 



tion 4.1 with noise cleaning applied. The distance up to 
which we will detect > 3 neutrinos with a probability of 



Elliptical E-SO < 0.05 

Spiral-like SOa-Sb 0.89 ± 0.33 

Spiral Sbc-Sd 1.68 ±0.60 

Others Sm, Irr., Pec. 1.46 =b 0.71 

Table 4: The expected rates of core-collapse supernovae 
for different galaxy types in supernova units (1 SNu = 
lSN(100yr)-i(10iOL|)-i). Values from ^, scaled by 
1.68 (see text). 



> 50% is found to be 6.9 Mpc for the TBP model, 9.8 Mpc 
for the LL model and 25.6 Mpc for the dark SN model, 
respectively. 

With this probability at hand along with the super- 
nova rate in the local environment, we can compute the 
expected supernova detections. We start from a catalog of 
nearby galaxies that reaches up to 100 Mpc [53]. Follow- 
ing [54] we assume that the blue luminosity of a galaxy is 
proportional to the star formation rate and hence also to 
the supernova rate. The conversion factor will depend on 
the galaxy type and is obtained from SN observations (cf. 
Table [4| . These conversion factors lead to a total SN rate 
that is lower by a factor 1.68 compared to a more recent 
result derived from a comprehensive compilation of local 
SNe [55 . We include this additional scale factor in our 
rate estimate to ensure consistency with currently avail- 
able data. For a limited number of galaxies the catalog 
leaves the type unspecified. Testing the allowed range of 
conversion factors, this leads to an error < 4% on the total 
core-collapse SN rate for distances up to 20 Mpc. 

Comparing with a theoretical prediction based on the 
initial mass function and cosmological star formation rate [9] , 
our SN rate is still a factor two lower. One explanation 
might be that the conversion factors in Table [4] are based 
on observations with a bias to miss many faint SNe. An- 
other possibility is a significant contribution of dark su- 
pernovae, that would not be detected optically, but still 
emit neutrinos [5. Scaling our blue luminosity prediction 
by a factor of two, we also find good agreement with the 
observed rate of nearby SNe [56] , and regard the rate esti- 
mates based on actual observations of SNe and those scaled 
to the star formation rate — both shown in Figure [TT] — as 
lower and upper limit, respectively. 

Table [5] gives a summary of expected SN detections per 
decade in different neutrino event multiplicity bins. The 
total resulting number of SN detections for the simulated 
127-string detector ranges from 23 to 46 per decade for 
the LL model and about half of that for the TBP model. 
For dark SNe we make the modest assumption that SNe 
collapse to a black hole at a rate of 10% of the regular 
core-collapse SN rate. Yet due to a more energetic neu- 
trino spectrum, dark SNe are detected at a higher effi- 
ciency resulting in between 20 and 40 dark observations 
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Figure 11: Range of cumulative expected supernova rate: 
lower curve (dashed) according to observed SN rate and 
blue luminosity [54], normalized to the results of [55], up- 
per curve (solid) scaled by a factor 2 to match expectations 
from the cosmic star formation rate [9] . The galaxy distri- 
bution and their blue luminosity values are taken from [53] . 



per decade, comparable to the number of detections from 
the LL model. Altogether, one can expect to observe at 
least one SN per year on average, perhaps up to 5 or more. 
The rate of SNe producing strong bursts of ten or more 
neutrinos is between 0.3 and 0.6 per year without dark 
SNe, reaching up to ^ 1 event per year if dark SNe are in- 
cluded. Remember that these numbers are obtained after 
applying the cuts to reject triggers from a sensor self- noise 
rate of 10 Hz as outlined in section |4]T| Cuts to fully reject 
the physical background of solar neutrinos are not yet in- 
cluded but we expect only a small reduction of rates from 
these. 

6. Conclusion 

A (9(10) megaton scale neutrino detector is required 
to extend the sensitivity to core-collapse SNe beyond the 
Large Magellanic Cloud to neighboring galaxies and yield 
a detection rate of up to several SNe per year. In this paper 
we have explored two implementations for such a detector 
in the Antarctic ice: either in the diffusive ice (around 
900 m below surface) where photon transport is diffusion 
dominated or in very clear ice (at around 2300 m depth) 
where photons can travel with very little scattering. In 
both cases, the detector geometry would resemble that of 
IceCube, however, with a much reduced string spacing and 
significantly increased density of photo-sensors per string. 
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CCSNe 


CCSNe 


Dark SNe 


neutrinos 


(LL) 


(TBP) 




3 




8.7 


3.6 


5.9 


4-5 




5.6 


2.8 


6.5 


6-9 




3.6 


2.3 


4.5 


^10 




5.5 


2.6 


3.4 



sum ^ 3 (L|) 
sum ^ 3 (SFR) 



23.4 
46.7 



11.3 
22.7 



20.3 
40.5 



Table 5: Expected number of supernova detections within 
one decade by observed neutrino multiplicity based on per- 
galaxy rates from blue luminosity [54^ 55 . The last lines 
show the sum of SN detections for this rate prediction 
(Lq), as well as for a prediction scaled to match the star 



formation rate (SFR) [9 (c.f. Figure 11). Dark SNe are 



assumed to occur at a fraction of 10% of all CCSNe. 



To achieve the (9(10) megaton scale with sensitivity to 10 
MeV neutrinos at trigger level would require a 61-string 
installation with about 20,000 photo-sensors similar to the 
HQE sensors employed by IceCube for the diffusive ice 
location and about ten times that total photo-cathode area 
for the clear ice location. 

Such a detector could be built in a similar manner as 
IceCube, by drilling holes into the ice and deploying strings 
holding the light sensors. However, with the demandingly 
large number of photo-sensors required, dedicated R&D 
for photo-sensors with large effective photo-cathode area 
seems appropriate. Besides the large effective area, we 
identify a low noise rate as crucial requirement for such 
a detector. Even for a self-noise rate of 10 Hz per meter 
instrumented string, the spatial and temporal distribution 
of hits has to be exploited and the number of strings in- 
creased to 127 to retain the (9(10) megaton target effective 
mass. The short scattering length will render this method 
void for the diffusive ice, making the task to separate neu- 
trino events from random noise triggers yet more difficult. 

The main physical background remaining after the self- 
noise reduction arises due to solar ^B neutrinos. Given a 
reconstruction algorithm, these can be identified via their 
angular proximity to the sun and lower photon multiplic- 
ity. Atmospheric muons will provide exceedingly bright 
events in such a dense installation, resulting in a small 
downtime for the deep ice. More challenging is the rejec- 
tion of spallation events and Michel electrons, yet with the 
IceCube detector fully surrounding the array, their iden- 
tification will be straight-forward and only a modest per- 
formance reduction will arise from their rejection. 

Using a catalog of nearby galaxies, we have computed 
the rate of detectable SN neutrino bursts. We show that 
depending on the SN rate and explosion model assumed, 
we can expect to observe between 11 and 46 SNe per 
decade in neutrinos (not counting the dark SNe) in the 
127 string detector configuration. By combining the ob- 
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servation of a neutrino burst with an optical detection of a 
SN, one can thereby disentangle the question of explosion 
model and SN rate. We note that in the future, with large, 
wide field optical surveys covering essentially the full sky, 
nearby SNe should be found in an even more systematic 
manner than today. 

Dark supernovae, where the star collapses to a black 
hole, could produce more than 20 detectable bursts per 
decade. Such a SN can be indirectly identified through the 
absence of an optical counterpart (with the risk of confus- 
ing it with a regular, dust obscured SN), or more directly, 
by observing neutrinos of higher energies. The limited en- 
ergy resolution of the detector (^ 30% per neutrino event) 
should be sufficient to estimate the effective temperature 
of the neutrino emission and hence the origin of the burst. 

A neutrino detector as described in this paper will not 
only yield a precise measurement of the local supernova 
rate and can uncover dark supernovae. A few of the su- 
pernovae will be closer, or perhaps even galactic, and hence 
yield a much higher number of coincident neutrinos, allow- 
ing to infer details about the explosion or even the neu- 
trino hierarchy. Lastly we point out that the low energy 
threshold and very large effective mass make such a detec- 
tor potentially interesting for a number of other physics 
phenomena, such as e.g. proton decay studies. While we 
have demonstrated that one can achieve the desired goal - 
routine observation of SNe in neutrinos - we acknowledge 
the assumptions we have made on our way. In particular, 
we have identified the development of large-area, low- noise 
photo-sensors suitable for deployment in the ice at South 
Pole as a key requirement. 
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